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The Multi-periodic Blazhko Modulation of CZ Lacertae' 
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ABSTRACT 

Thorough analysis of the multicolour CCD observations of the RRab-typc variable, 
CZ Lacertae is presented. The observations were carried out in two consecutive ob- 
serving seasons in 2004 and 2005 within the framework of the Konkoly Blazhko Survey 
of bright, northern, short-period RRab variables. The O — C variation of CZ Lac in- 
dicated that a significant period decrease took place just around the time of the CCD 
observations. Our data gave a unique opportunity to study the related changes in the 
pulsation and modulation properties of a Blazhko star in detail. Two different pe- 
riod components (« 14.6 d and ps 18.6 d) of the Blazhko modulation were identified. 
Both modulation components had similar strength. The periods and amplitudes of 
the modulations changed significantly from the first season to the next, meanwhile, 
the mean pulsation amplitude slightly decreased. The modulation frequencies were 
in a 5:4 resonance ratio in the first observing season then the frequencies shifted in 
opposite directions, and their ratio was close to the 4:3 resonance in the next season. 
The interaction of the two modulations caused beating with a period of 74 d in the 
first season, which resembled the 4-yr-long cycle of the sa 40-d modulation of RR Lyr. 
The mean values of the global physical parameters and their changes with Blazhko 
phase of both modulation components were determined by the Inverse Photometric 
Method. 

Key words: stars: horizontal branch - stars: variables: other - stars: individual: 
CZ Lac - stars: oscillations (including pulsations) - methods: data analysis - tech- 
niques: photometric 



1 INTRODUCTION 

RR Lyrae variables have great astrophysical importance for 
they are tracers of stellar evolution, they are used as distance 
indicators and they are test objects of stellar pulsation theo- 
ries. We think that these objects are known fairly well. They 
are pure radial pulsators oscillating in fundamental mode, 
in first overtone and, in some cases, in both. 

There is, howeve r, a flaw in th is nice picture: the so- 
called Blazhko effect (|BIazhkdll907h , the periodic variation 
of the light-curve shape of some of the RR Lyrae stars. Al- 



* Based on observations collected with the automatic 60-cm tele- 
scope of Konkoly Observatory, Budapest, Svabhegy 
f E-mail: sodor@konkoly.hu 



though the phenomenon has been known for more than a 
century, its physical origin is still unknown. Recently, it 
has been found that a much larger percentage of the RR 
Lyrae stars exhibi t light-curve modula t ion than it was pre - 
viously assumed l|jurcsik et al.l [20 09b: Bc nko et alj|201(ih . 
which makes the problem even more embarrassing. 

To explain the light-curve modulation, several mod- 
els have been put forward during the past decades, but 
none of them is able to reproduce most of the ob- 
served properties of the modulation. The oblique mag- 
netic rotator models (|Cousensl 1 19851 : IShibah ashi 200j|) are 
based on the presumption of a bipolar magnetic field 
with a magnetic axis oblique to the axis of rotation . 
The resonance models dNowakowski fc Dziembowskil 1200 ll : 
iDziembowski fc Mizerskl 2004 ) explain the modulation with 
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Table 1. Log of the CCD observations of CZ Lac obtained with 
the 60-cm automatic telescope. Note that, due to technical prob- 
lems, observations were not obtained in each band on 10 nights 
in the second observing season. 



Season 


From [JD] 


To [JD] 


Filter 


Nights 


Points 


1st 


2453 266 


2 453 411 


B 


68 


3669 








V 


68 


3683 








Rc 


68 


3649 








Ic 


68 


3629 


2nd 


2 453 648 


2 453 731 


B 


12 


3011 








V 


52 


4460 








Rc 


50 


3446 








Ic 


50 


3371 



the excitation of non-radial modes o f low radial deg r ee. Ac - 
cording to the model proposed by IStothersI (|2006l . l2010t ). 
a turbulent dynamo in the envelope triggers the light-curve 
modulation. The analysis of the 127-d-long quasi-continuous 
Kepler data raised the possibility that a 9:2 radial-overtone 
- fundamental-mode r esonance might b e connected to the 
Blazhko phenomenon ()Szab6 et al.ll201ol ). 

To find the physical explanation for the Blazhko effect, 
we have to know how it manifests itself in the observed 
properties, i.e. detailed observations of Blazhko stars are 
needed. Therefore, w e launched the Konko ly Blazhko Survey 
in 2004 jSodorl 120071 ; lJurcsik et all 1200901 ). which aimed at 
collecting extended, accurate, multicolour photometric data 
of Blazhko stars. In this paper, we report the results on 
CZ Lacertae, one of the most intriguing target of our sur- 
vey. 

CZ Lac (m v » ll m 7, Q2000 = 22 h 19 m 30 i !76, <5 20 oo = 
+51°28'14'.'8) is a bright, short-period, fundamental-mode 
RR Lyr ae-type (RRa b) variable in the northern sky discov- 
ered bv lFlorial (|l940h . No photometric study of CZ Lac has 
been published yet. 



2 THE DATA 
2.1 Observations 

CZ Lac has an optical companion at 10-arcsec separation, 
which makes low-resolution photometric observations defec- 
tive. 

CZ Lac was observed with a photometer attached to 
the 60-cm telescope of the Konkoly Observatory, Svabhegy, 
Budapest, through B and V filters on seven nights in 1967. 
These observations showed light-curve variability, however, 
no secondary period were found. The light-curve changes 
were attributed to the light contamination of the close com- 
panion at that time. The photoelectric data are available 
online as Supporting Information (Table SI). This table lists 
the differential B and V magnitudes of CZ Lac relative to 
the comparison star, BD +50°3664. 

We observed CZ Lac with the same, refurbished and 
automated 60-cm telescope equipped with a Wright Instru- 
ments 750 x 1100 CCD camera (field of view 17 x 24 ar- 
cmin) in 2004 and 2005. About 29 000 frames were obtained 
in BV(RI)c bands. The data of the two seasons spanned 
146 and 84 d. The gap between the consecutive seasons' ob- 
servations was 237-d-long. Due to problems with the filter 



Table 2. The standard reference magnitudes of CZ Lac corrected 
for atmospheric extinction in the reference images and their esti- 
mated uncertainties are given in the table. 



Filter 


magnitude 


B 


11.83 ±0.02 


V 


11.50 ±0.01 


Rc 


11.25 ±0.01 


Ic 


10.96 ±0.03 



changer, B, Rc and Ic observations were not obtained on 
each of the nights in the second season. The log of the CCD 
observations is presented in Table Q] 

To tie the instrumental CCD magnitudes to the stan- 
dard Johnson-Cousins system, observations of stars in a 
20 x 20 arcmin field centred on CZ Lac were made with the 
USNO Flagstaff Station 1.0-m telescope equipped with a 
SITe/Tektronix 1024 x 1024 CCD. The positions and the 
BV(RI)c magnitudes of these stars are available online as 
Supporting Information (Table S2). 

2.2 CCD reduction, photometry 

Standard CCD calibration was done using the program 
package IRAlfl In order to eliminate the light contamination 
of the close compa nion, Image S ubtraction Method (ISM), 
as implemented by lAlardl (|2000t) , was applied on the cali- 
brated frames. 

The relative fluxes of the variable, i.e. the differences 
between the measured fluxes and the flux of the variable 
in a reference image, were determined for each frame by 
aperture photometry using the daophot package of iraf. 
To convert the relative fluxes to differential magnitudes, 
the fluxes and the corresponding magnitudes of the vari- 
able were measured in the reference image by aperture 
photometry using the standard magnitudes of 3 surround- 
ing field stars (USNO-A2.0 1350-16354156, ...16333516 and 
...16335838, see Table S2). The reference magnitudes cor- 
rected for atmospheric extinction and transformed into the 
standard Johnson-Cousins system are listed in Table [2] The 
standard relative BV(RI)c time series of CZ Lac are given 
in electronic form as Supporting Information (Tables S3- 
S6). The differential magnitudes correspond to the reference 
magnitudes given in Table [2] 



3 PULSATION-PERIOD VARIATIONS: THE 
O-C DIAGRAM 

Normal maximum times were derived for the middle of the 
observations of the two seasons' CCD V band data by fitting 
mean pulsation light curves. These normal maximum times 
and the times of the 6 photoelectric V maxima observed in 
1967 are listed in Tabled 

There are 87 light maxima of CZ Lac collected in the 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 3. Maximum times of the Konkoly photoelectric and CCD 
observations. 



HJD 


observation 


2439740.437 


photoelectric 


2439754.272 


photoelectric 


2439758.587 


photoelectric 


2439762.478 


photoelectric 


2439765.502 


photoelectric 


2439766.360 


photoelectric 


2453338.647 


CCD 2004 normal max. 


2453689.567 


CCD 2005 normal max. 



data base (see references therein). We omitted 4 dis- 
crepant, visual observations and complemented the data set 
with the 8 Konkoly data. The O — C diagram, constructed 
from these maximum times, is plotted in the top panel of 
Fig II The following ephemerida were used for the calcula- 
tion: 

t£XL c) [HJD] = 2 453 338.647 + d 432173 E. 

The initial epoch and the period were taken from the light- 
curve solution of the first season's CCD observations. There 
might be uncertainties in the epoch counts when large gaps 
occurred between the observations, for example, around 
HJD 2 447 000, when observations were missing for a 3300- 
d-long interval. 

The O — C diagram shows that the pulsation period of 
CZ Lac underwent a significant decrease just about when 
the Konkoly CCD observations were obtained. The O — C 
plot after 1990 is shown magnified in the middle panel of 
Fig [2 These data are fitted linearly before and after 2005.5 
(HJD 2 453 520); the corresponding pulsation periods are 
0.432183 and 0.432138 d. The pulsation-period changes can 
be followed in the bottom panel of Fig. [TJ Periods deter- 
mined from the O — C data prior to and after 2005.5 and 
found in the two seasons' CCD observations are shown. Our 
CCD data indicates that a small period decrease occurred 
between the two observing runs (see Table[4|. As the periods 
determined for 2004-2005 fell between the two period values 
obtained from the O — C data, the pulsation-period change 
of CZ Lac should not have bee have abrupt. The period- 
change event between the extrema should have lasted for 
several years. 

The pulsation period of Bl azhko stars ofte n shows 
abrupt, ir r egular changes (|Le Borgne et all 120071 ; 
ISzeidl et all |2010|) . Evidences have also been given 
that changes in the modulation properties ar e connected 
to th e pulsation-period variations (e.g., in ISodor et al.l 
2007). However, no continuous, extended observation of 
any Blazhko star have ever been obtained coincidentally 
with a rapid period change event. CZ Lac was observed in 
2004-2005, just around the time when a rapid change in 
its period began. It gave us a serendipitous opportunity to 
study such an event in detail. 
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Figure 1. O - C diagram of CZ Lac between 1937 and 2009. 
Data from the literature (collected in the CEOS data base) are 
plotted with dots, Konkoly observations are plotted with dia- 
monds. The data after HJD 2 448 000 (1990) are shown magnified 
in the middle panel. These data indicate a sudden period decrease 
just around the time of the CCD observations (2004-2005). The 
straight lines fitted to the data shown in this panel prior to and 
after 2005.5 (HJD 2 453 520) are also drawn. Bottom panel shows 
the period values determined from the O — C data and from the 
Konkoly CCD observations. This figure proves that, contrary to 
the O — C results, the pulsation period of CZ Lac did not change 
abruptly, but the period change lasted for several years. 



4 LIGHT-CURVE ANALYSIS 
4.1 Light-curve solutions 

The CCD V light curves of the two observing seasons, folded 
with the pulsation period, are shown in Fig. [2] The plots in- 
dicate that the modulation properties changed notably be- 
tween the two seasons. The overall strength of the amplitude 
and phase modulations decreased, which did not allow us to 
analyse all the data together. Therefore, the two seasons 
were treated separately. 

A mathematical description of the brightness varia- 
tions was sought in the form of Fourier sums, utilising 
discrete Fourier transformation and non-linear and linear 
least-squa res fitting met hods. We used the program package 
mufran (|Kollath|[l990l '). 



the fitting abilities of GNUPLOl 



2 http://rr-lyr.ast.obs-mip.fr/dbrr/dbrr-Vl.0J3.php 3 http://gnuplot.info/ 
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Figure 2. CCD V light curves of CZ Lac folded with the pulsation period for the two observing seasons. The properties of the modulation 
changed significantly between the two seasons. The overall strength of both the amplitude and phase modulations decreased. 



and an own-developed code, which fitted the data with in- 
dependent frequencies and their linear combinations non- 
linearljjf). 

The Fourier components were identified in the data suc- 
cessively until the residual spectra showed peaks higher than 
4<r in at least two colour bands at the same frequency. In 
addition, peaks at linear-combination frequencies of inde- 
pendent frequencies found in the spectra were taken into 
account down to about the 2a significance level, if they ap- 
peared in at least two bands of the same season's data. 

Different components at different frequencies were 
needed to accurately fit the two seasons' data, but the same 
components were used for the B, V , Rc and Ic data sets 
in one season. The solutions of the first and second seasons' 
data comprised 134 and 119 frequencies, respectively. 

The analysis revealed modulation components located 
symmetrically around the pulsation frequency (/ p ) and its 
harmonics in the Fourier amplitude spectra, which is char- 
acteristic of Blazhko stars. Side frequencies corresponding to 
modulations of two different frequencies (/ m i and / m 2) were 
detected in both seasons. All the identified Fourier compo- 
nents were derived from the three independent (base) fre- 
quencies: / P ,/mi and / m 2- Note that the 'ml' and 'm2' in- 
dices are arbitrary, and do not reflect the dominance of the 
/ m i modulation component. 

Several modulation peaks appeared in the spectra 
around the pulsation peaks with separations of the two mod- 
ulation frequencies and their linear combinations. The mod- 
ulation frequencies f m i and / m 2 were also identified in the 
spectra. The most intriguing feature is the detection of the 
frequencies with separations of half of / m 2 (kf p ± 0.5/ m 2). 
Linear combination terms involving this half frequency also 
appeared in the first seasons' data. 

Having the relevant Fourier components identified, their 
frequency values were fixed to the values obtained from the 
V^-band light curve solutions as the signal-to-noise properties 
were the best in this band. The V-band solutions were ob- 
tained by non-linear fits that allowed to adjust the three in- 
dependent base frequencies, and the combination frequencies 
were locked to the base-frequency values during the fitting 
process. The base frequencies, their absolute and relative 
changes and the formal la errors of the V-band solutions 
are given in Table [4] for the two seasons. 



4 http : / /konkoly. hu /staff /sodor /nlfit . html 



To determine the amplitudes and phases of the fre- 
quency components in the four bands, linear fits were calcu- 
lated with all the frequencies fixed to their accepted values. 
The parameters of the solutions are listed in Tables [5] and [5] 
The formal la errors of the Fourier amplitudes and phases 
are also given. 

The detected frequencies are illustrated schematically 
in Fig. [3] The complexity of the multiplet structures in the 
Fourier spectra, and the changes in the frequencies, ampli- 
tudes and in the appearing linear combination terms be- 
tween the two seasons can be well followed in this figure. 

The residual spectra showed no signal stronger than 
4cr, except one peak in each of the B, Rc and Ic bands, 
which appeared, however, at different frequencies in the 
first season's data. These peaks were most probably arte- 
facts and were not intrinsic to the variable. Inspecting the 
residual light curves folded with the pulsation period, we 
found extended scatter around the phase of the middle of 
the rising branch. Similar beh aviour of the r e siduals of other 
Blazhko stars were shown in lJurcsik et all (|2008al ) and in 
IPoretti et ah! ||201Ch . 

The Fourier frequency components kf p ± 0.5/ m 2 are 
extraordinary. No similar peak has ever been detected in 
any other Blazhko star. A recent study on the long-term 
behaviour of Blazh ko stars in the globular cluster, M5 
(|jurcsik et al.l 120101 1 revealed that V2/M5 had a modula- 
tion period of 67 d at around the middle of the 20th cen- 
tury, while its Blazhko cycle turned to be 130 d, just about 
twice of its previous value, after 1960. Anomalous modu- 
lation components with 11.5,f m , 12.5.f m and 13.5/ m separa - 
tions were also detected in MW Lyr |jurcsik et al. 2008a). 
These examples indicate that the half of the modulation fre- 
quency might have some relevance in the interpretation of 
the data. It cannot even be excluded that 0.5/ m 2 is one of 
the real modulation frequencies of CZ Lac and peaks with 
separations of / m 2, 2/ m 2 and 3/ m 2 are, in fact, its 2nd, 4th 
and 6th modulation harmonics. The side-lobe peaks with 
0.5/ m 2 separations around the pulsation harmonics are, how- 
ever, much weaker than those with / m 2 separations. There- 
fore, keeping in mind the possibility that the true modu- 
lation frequency might be 0.5/ m 2, for practical purposes we 
have chosen the large-amplitude component, / m 2, as the first 
harmonic-order component of this modulation. 
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Table 4. Independent (base) frequencies of the light-curve solutions (/ p - pulsation frequency; / m i, / m a - modulation frequencies), their 
errors and their absolute and relative changes between the two seasons. The last two columns list the corresponding period values. 



1st season 2nd season Preq. change 1st season 2nd season 

component frcq. [cd^ 1 ] freq. [cd _1 ] absolute [cd — 1 ] relative Period [d] Period [d] 



/p 


2.313887 


2.313910 


2.3 


10" 


-5 


9.9 


io- 


6 


0.4321732 


0.4321689 


± 


0.000002 


0.000003 


0.5 


10" 


-5 


2.2 


io- 


6 


0.0000004 


0.0000006 


/ml 


0.05406 


0.05347 


-5.9 


io- 


-4 


-1.1 


10" 


2 


18.50 


18.70 


± 


0.00002 


0.00007 


0.9 


io- 


-4 


0.2 


10" 


2 


0.01 


0.02 


/m2 


0.06757 


0.06954 


19.7 


io- 


-4 


2.9 


10" 


2 


14.804 


14.38 


± 


0.00002 


0.00005 


0.7 


io- 


-4 


0.1 


10" 


2 


0.004 


0.01 



Table 5. Frequencies, amplitudes (in mag) and phases (in rad) of the B, V, Rc and Iq light-curve solutions of CZ Lac in the first 
observing season. The full table is available as Supporting Information with the online version of this article. 







B 




V 




Rc 




Ic 




identification 


/ [ccT 1 ] 


amplitude 


phase 


amplitude 


phase 


amplitude 


phase 


amplitude 


phase 


1/p 


2.313887 


0.4894(6) 


3.393(1) 


0.3535(5) 


3.330(1) 


0.2763(4) 


3.248(2) 


0.2110(5) 


3.117(2) 


2/p 


4.627774 


0.2537(6) 


2.909(2) 


0.1895(5) 


2.896(2) 


0.1500(5) 


2.868(3) 


0.1137(5) 


2.822(4) 


3/p 


6.941660 


0.1472(6) 


2.737(4) 


0.1115(5) 


2.731(4) 


0.0891(5) 


2.726(5) 


0.0691(6) 


2.728(7) 


/ml 


0.054065 


0.0056(4) 


1.51(7) 


0.0035(3) 


1.40(10) 


0.0035(3) 


1.50(9) 


0.0035(3) 


1.84(9) 


/ma/2 


0.033786 


0.0025(4) 


0.8(2) 


0.0024(3) 


0.5(1) 


0.0029(3) 


0.9(1) 


0.0040(3) 


0.87(9) 


/m2 


0.067573 


0.0041(4) 


3.0(1) 


0.0013(3) 


2.5(2) 


0.0027(3) 


2-3(1) 


0.0023(4) 


1.4(1) 


1/p — /ml 


2.259822 


0.0235(5) 


1.34(2) 


0.0172(4) 


1.33(2) 


0.0137(4) 


1.41(3) 


0.0107(4) 


1.52(4) 


2/p — /ml 


4.573709 


0.0272(5) 


0.67(2) 


0.0202(4) 


0.69(2) 


0.0159(4) 


0.74(2) 


0.0122(4) 


0.78(3) 


1/p + /ml 


2.367951 


0.0515(5) 


3.36(1) 


0.0378(4) 


3.41(1) 


0.0289(4) 


3.46(1) 


0.0223(4) 


3.53(2) 


2/p + /ml 


4.681838 


0.0448(5) 


2.73(1) 


0.0331(4) 


2.77(1) 


0.0270(4) 


2.78(2) 


0.0208(5) 


2.79(2) 


1/p - /m2 


2.246314 


0.0136(4) 


0.85(3) 


0.0105(3) 


0.84(3) 


0.0091(3) 


0.83(4) 


0.0071(4) 


0.79(5) 


2/p - /m2 


4.560201 


0.0160(4) 


0.25(3) 


0.0117(3) 


0.33(3) 


0.0102(3) 


0.37(3) 


0.0080(4) 


0.35(4) 


1/p + /ma 


2.381459 


0.0391(4) 


4.99(1) 


0.0291(3) 


5.03(1) 


0.0229(3) 


5.00(1) 


0.0165(4) 


4.96(2) 


2/p + /m2 


4.695346 


0.0316(4) 


4.37(1) 


0.0238(4) 


4.39(1) 


0.0179(3) 


4.42(2) 


0.0144(4) 


4.54(3) 



Table 6. Frequencies, amplitudes (in mag) and phases (in rad) of the B, V, Rc and Ic light-curve solutions of CZ Lac in the second 
observing season. The full table is available as Supporting Information with the online version of this article. 



identification 


/ [cd- 1 ] 


B 

amplitude 


phase 


V 

amplitude 


phase 


Rc 

amplitude 


phase 


Ic 

amplitude 


phase 


1/p 
2/p 
3/p 


2.313910 
4.627821 
6.941731 


0.4854(5) 
0.2495(5) 
0.1454(5) 


2.863(1) 
1.843(2) 
1.131(4) 


0.3509(4) 
0.1863(4) 
0.1109(4) 


2.802(1) 
1.830(2) 
1.121(3) 


0.2749(4) 
0.1469(4) 
0.0886(4) 


2.721(1) 
1.808(3) 
1.132(4) 


0.2046(4) 
0.1094(4) 
0.0689(4) 


2.577(2) 
1.783(4) 
1.139(6) 


/ml 

/m2 
1/p — /ml 
2/p — /ml 


0.053469 
0.069541 
2.260441 
4.574351 


0.0011(4) 
0.0036(4) 
0.0088(6) 
0.0094(5) 


4.0(4) 
2.9(1) 
4.17(8) 
3.15(6) 


0.0015(3) 
0.0018(3) 
0.0055(4) 
0.0062(4) 


4.1(2) 
3.7(2) 
4.59(9) 
3.25(7) 


0.0015(3) 
0.0018(3) 
0.0039(4) 
0.0041(4) 


4.0(2) 
4.7(2) 
4.3(1) 
3.2(1) 


0.0014(3) 
0.0016(3) 
0.0033(5) 
0.0032(4) 


4.2(2) 
4.4(2) 
4.3(2) 
3.5(1) 


1/p + /ml 
2/p + /ml 


2.367380 
4.681290 


0.0191(6) 
0.0137(5) 


5.33(3) 
4.26(4) 


0.0151(4) 
0.0120(4) 


5.19(3) 
4.18(3) 


0.0118(4) 
0.0107(4) 


5.24(4) 
4.27(4) 


0.0094(5) 
0.0082(5) 


5.42(5) 
4.41(6) 


1/p — /m2 
2/p - /m2 


2.244370 
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Figure 3. Multiplet structures of the light-curve solutions of CZ Lac in the two observing seasons. Horizontal axes show the pulsation 
harmonic order (fc). The vertical axes show the vicinity of the pulsation harmonic frequencies. The identification of the frequency com- 
ponents are given. The V amplitudes arc colour-coded; amplitudes higher than 0.03 mag arc shown full black. These figures demonstrate 
that significantly different frequency components are needed to fit the two seasons' data. Note the discernible changes in the amplitudes 
and separations of the modulation components. 



4.2 Changes between the seasons 

It was shown in the previous sections that a slight increase of 
the pulsation period (0.0000043±0.0000010 d) was accompa- 
nied by significant frequency and amplitude changes of the 
complex modulation of CZ Lac. The CCD V magnitudes 
of the two seasons are plotted against HJD in Fig. [4] The 
upper envelopes derived from the synthetic light curves are 
also drawn. The interaction between the two modulations 
produced beating in the first season; they cancelled each 
other out periodically as can be seen around HJD 2 453 290 
and HJD 2 453 365. The envelope of the second season shows 
a weaker amplitude modulation. 

The kf p ± / m i modulation components had somewhat 
higher amplitudes than the amplitudes of the kf p ± / m 2 fre- 
quencies in the first season. The amplitudes of both modu- 
lations decreased to the second season by different rates, so 
that the amplitude relation of the two modulations reversed. 
The amplitudes of the kf p ± / m 2 modulation components 
were larger than those of the kf p ± / m i components then. 

It seemed that no changes occurred in the frequencies 
and amplitudes within one season, so the pulsation and 
modulations of CZ Lac could be described reasonably well 
by stationary harmonic functions in both seasons. The first 
season's data were extended enough to be divided into two 
parts at HJD 2 453 350. A separate analysis indicated that no 
systematic changes in the frequencies and / or amplitudes oc- 
curred in the first observing season. Had the amplitudes and 
frequencies changed somewhat during any of the observing 



seasons, the residual spectra would show significant peaks 
around the subtracted frequencies, which was not the case. 
These arguments support that the changes took place dom- 
inantly between the two seasons, within a relatively short 
interval of 237 d. 

Not only the amplitudes of the modulation components 
but also the amplitude of the pulsation changed between 
the two seasons. The peak-to-peak mean pulsation ampli- 
tudes of both seasons were determined from the synthetic 
light curves taking into account only the pulsation harmon- 
ics. Their errors were estimated as the sum of the errors of 
the amplitudes of the pulsation harmonic components. Be- 
cause the errors of the amplitudes were correlated and the 
phases were not taken into account, these errors overesti- 
mated the true uncertainties. We found that the mean pul- 
sation amplitude decreased by 0.027 ± 0.011, 0.023 ± 0.009, 
0.015 ± 0.008 and 0.017 ± 0.009 mag in B, V, R c and 7 C 
bands, respectively, simultaneously with the diminution of 
the modulation amplitudes. The epoch-independent phase 
differences of the pulsation harmonics did not change signif- 
icantly between the two seasons, indicating that the shape 
of the mean pulsation light curve remained stable, whilst its 
amplitude decreased. 

Simultaneous pulsation-period and -amplitude changes 
were also detecte d in RR Gem around the year 1938 
l|S6dor et all 120071 ), when the character of the modulation 
of the star changed as well. The pulsation period and ampli- 
tude increased in RR Gem, while they decreased in CZ Lac. 
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Figure 4. CCD V light curves of CZ Lac versus HJD. The upper envelope of the data shows the brightness variations of the light maxima. 
The properties of the modulation changed considerably between the two seasons. The overall strength of the modulation decreased. 



4.3 Resonances; one or two modulation 
frequencies? 

The frequency ratio of the two modulations was / m 2//mi = 
1.2499 ± 0.0008 in the first season, i.e. the frequencies were 
in a 5:4 resonance ratio well within the error limit. As Ta- 
bleUshows, the modulation frequencies changed in opposite 
directions between the two seasons, and as a consequence es- 
caped from the 5:4 resonance. The frequency ratio of the two 
modulations was 1.301 ± 0.002 in the second season, which 
was closer to the 4:3 resonance ratio than to the 5:4 one. 

The transition from the 5:4 resonance to the near 4:3 
resonance hints that the two modulation frequencies prefer 
to be in resonance. Considering this finding, one might won- 
der if indeed two linearly independent modulation frequen- 
cies exist in CZ Lac. Could the light curves be described 
with only one modulation frequency and many of its har - 
monics, as, for ex ample, in MW Lyr |jurcsik et all [2008al ) 
and in V1127 Aql (jChadid et al.ll2010h ? In this case, the base 
modulation frequencies of the two seasons were significantly 
different, they were / m = 0.006757 cd" 1 = / m i/8 = / m2 /10 
and / m = 0.008802 cd" 1 « / m i/6 w / m2 /8 in the first and 
second seasons, respectively. The residual scatters of the one- 
modulation-frequency fits are the same in the first season, 
and they are only 3-5 per cent larger in the second season 
than those of the two-modulation-frequency solutions. The 
one-modulation solution would also give a natural explana- 
tion for the appearance of the combination terms. However, 
serious problems emerge when one tries to interpret this so- 
lution. 

One of the problems is that it is ambiguous how to 
relate the different modulation components of the two sea- 
sons' solutions in this case. The one-modulation-frequency 
solutions incorporate different-order modulation sidelobe 
frequencies (kf p ± n/ m ) in the two seasons. Components 
corresponding to n = {2, 3, 5, 8, 10, 16, 18, 20} and n = 
{2, 4, 6, 8, 12, 14, 16, 24} are involved in the first and sec- 
ond seasons' solutions, respectively. In contrast, if we accept 
the two-modulation-frequency solution, a consistent identi- 
fication of the modulation components can be set, and the 
connection between the frequencies of the two seasons is 
straightforward . 

Another problem of the one-modulation-frequency so- 
lution is that, in this case, the modulation harmonics of 



CZ Lac would behave rather differently than those in other 
Blazhko stars. In well-observed Blazhko stars with quin- 
tuplets and higher-order multiplet structures detected, the 
frequency components at the smallest separation (n = 1) 
had always the larg est amplit ude (for example, V1 127 Aql - 
IChadid et aT]|201oi, D M Cyg - lJurcsik et "al1l2009al MW Ly: 



ljurcsik et al.ll2008al . RV UMa - iHurta et al.ll2008h . In con- 
trast, according to the one-modulation solution of CZ Lac 
the largest-amplitude, dominant modulation peaks are the 
n — 8 and 10 components in the first season and the n = 8 
and 6 components in the second one. The n = 1 component 
is not ev en detectable in any o f the seasons. It has been 
shown in ISzeidl fc Jurcsikl <|2009h that the / p ± f m frequen- 
cies are always the largest-amplitude modulation compo- 
nents in the frequency spectrum of an amplitude- and phase- 
modulated harmonic oscillation, provided that the ampli- 
tude of the phase modulation is relatively small compared 
to the period of the oscillation. Note that the same argu- 
ment applies against the acceptance of 0.5/ m 2 as the base 
modulation frequency of the second modulation. 

All these arguments lead us to the conclusion that in- 
deed two independent modulations of CZ Lac coexisted and 
their frequencies changed in opposite directions between the 
two observing seasons. 



5 RESEMBLANCE BETWEEN THE 

MODULATIONS OF CZ LAC AND RR LYR 

The f» 40-d-long Blazhko mod ulation of RR Lyr s eemed 
to be ceased in every 4 years (jDetre fc Szei dlll973h. The 
latest weak-m odulation phase had been detected in 1975 
(jSzeidll Il976l ) in the Konkoly observations, which contin- 
ued until 1981. Although RR Lyr was not regularly ob- 
served in the past three decad es, photometric o bserva- 
tions were obtained in 1990- 1993 (ISzeidl et allll997T ), 1993- 
1994 dBelserenelll997ri. 1996 ijSmith et al]l2003h. 2003-2004 
dKolenberg et al.l l200rj ) and 2006-2007 (jKolenberg et all 
l2008h . All these observations showed maximum-light vari- 
ation in V band of at least 0.2-mag amplitude during 
the Blazhko cycle. The amplitude of the modulation was 
also large in 2009 accordin g to the recent Kepler ob- 
servations (|Szab6 et all l2010j). The spectrosco pic observa- 
tions from 1994-1995 (|Chadid fc Gilletl ll996T ), 1994-1997 
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Table 7. Mean physical parameters of CZ Lac and their errors derived from the mean light curves by the inverse photometric method, 
using static atmosphere models of [Fe/H] = —0.20 metallicity. 



M 


d 


E(B-V) 


M v 


(B - V) 


(V - I) 


R 


L 




[M ] 


[pc] 


[mag] 


[mag] 


[mag] 


[mag] 


[R@] 




[K] 


0.60 ±0.01 


1188 ±8 


0.19 ±0.03 


0.66 ± 0.01 


0.276 ± 0.002 


0.344 ± 0.002 


4.46 ±0.03 


46.5 ±0.5 


7100 ± 10 



1st season 2nd season 




0.5 1-0.5 0.5 0.5 1-0.5 0.5 




0.5 1 -0.5 0.5 0.5 1 -0.5 0.5 

Pulsation phase Blazhko phase Pulsation phase Blazhko phase 



Figure 5. Synthetic V light curves of CZ Lac corresponding only one of the modulation components are phased according to the 
pulsation and modulation frequencies for the two seasons. These data were used as input for the IPM to determine the variations of 
the pulsation-averaged physical parameters with Blazhko phase. Lines indicate the phase and brightness variations in maximum and 
minimum light. 



Chadidl l200Ch. 1996-1997 (|Chadid et all ll99Sh and 2000 



Chadid fc Chapellierl lioOrj ). all showed clear evidence of 



strong amplitude modulation. It seems that, in spite of 
the above mentioned observational activities, the weak- 
modulation states of RR Lyr have been either missed or the 
cyclic behaviour of the modulation has been less pronounced 
during the last decades. 

Variati on in the modulation period o f RR Lyr was also 
reported by iKolenberg et alj (|200fj . l200l ). In contrast with 
the 40-41 d Blazhko period found in previous investigations, 
a 38-39 d period was detected in these recent studies. 

Many aspects of the Blazhko behaviour of CZ Lac re- 
sembles the properties of RR Lyr. The modulation of CZ Lac 
seemed to cease in every 74 days in the first season, as a 
consequence of the interaction between two, similarly strong 
modulations. This beating is shown by the envelope in Fig. [4] 
The beating period is l/(/ m 2 — /mi) = 74.02 d. We suppose 
that the four-year cycle of RR Lyr can also be explained by 
the beating of two modulations of similar strength, if the 
periods of the modulations differ by about 1 d. To resolve 
modulation peaks being so close to each other, continuous 
observations of more than four years are needed, which will 
be available from the Kepler mission. 

The strength of the modulation components of CZ Lac 
changed significantly from the first observing season to the 



next; while the 18.5-d modulation (/ m i) had somewhat 
larger amplitude than the 14.8-d modulation (/ m 2) in the 
first season, the amplitude of the longer-period component 
(/mi) was extremely small in the second one. Consequently, 
no cancellation (beating) of the modulations was detected 
then. A similar behaviour might explain the vanishing of the 
4-yr cycles of RR Lyr during the last decades and the strong 
change in its observed modulation period. 



6 CHANGES IN PHYSICAL PARAMETERS 

We have develo ped an inverse photometric Ba ade-Wesselink 
method (IPM) (|S6dor, Jurcsik fc Szeid1l2008l ) to determine 
the mean physical parameters of RR Lyrae stars, such as ef- 
fective temperature (T e g), absolute visual magnitude (My ), 
radius (R), distance (d), mass (M) and dereddened colours 
[(B — V)o, (V— ic)o]j from multicolour photometric observa- 
tions. It has been shown that using good quality multicolour 
photometry, IPM yields results of similar accuracy as di- 
rect spectroscopic Baade-Wesselink methods. The IPM also 
gives the variations with pulsation phase of those param- 
eters that change during the pulsation cycle. This method 
has been successfully applied to derive the mean physical pa- 
rameters and their changes during the Blazhko cycle already 
for four modulated RRab stars (MW Lyr - ljurcsik et alj 
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l2008bl, DM C yg - ljurcsik et all |2009al , RR Gem, SS Cnc 
- IS6doi|[2009l ). The IPM is a useful tool to unravel astro- 
physics buried deep in the hundreds of Fourier coefficients 
obtained by the light-curve-fitting process. 

In order to determine the changes with any of the mod- 
ulations of CZ Lac by the IPM, the two modulations had to 
be separated first. It was achieved by generating synthetic 
light curves from the solutions of each season's data taking 
into account only the pulsation frequencies and terms that 
belonged to one of the modulations only. The coupling mod- 
ulation terms were omitted, as they were not periodic with 
any of the modulations. Considering their small amplitudes, 
the omission of these terms presumably had a negligible ef- 
fect on the results. The synthetic V light curves of the two 
modulations for both seasons are plotted in Fig. \5\ 



Table 8. Differential magnitude zero-points of the solutions for 
the two observing seasons in the four photometric bands. Relative 
values with respect to the reference magnitudes listed in Table [2] 
are given. 



band 


1st season 


2nd season 


B 


0.2902(2) 


0.2907(3) 


V 


0.1538(3) 


0.1544(3) 


H c 


0.1327(3) 


0.1346(2) 


Ic 


0.1073(3) 


0.1086(2) 



same within the millimagnitude limit in the two seasons. 
The differential magnitude zero-points of the light-curve so- 
lutions for the two seasons are compared in Table [S] 



6.1 Constant parameters 



6.2 Blazhko-phase dependent parameters 



Following the methodology of our previous investigations, 
we started the analysis by determining those parameters 
that obviously did not change. These were the metallicity 
([Fe/H]), mass (SOT) , distance (d) and the interstellar red- 
dening [E(B - V)] of the star. 

T he metallicity of CZ Lac was determ ined spectroscop- 
ically (ISuntzeff. Kraft fc Kinmannl Il994f ) and photometri- 
cally (|jurcsik et al.ll2009bl ) to be —0.13 and —0.26, respec- 
tively. Therefore, static atmosphere model grids c orrespond- 
ing to [Fe/H] = -0.20 |Castelli fc Kurucj]2003T ) were used 
for the IPM analysis. 

With the exception of the metallicity, all the other 
constant parameters were derived by the IPM using the 
mean pulsation light curves, which were generated from the 
light-curve solutions taking into account only the pulsa- 
tion components. To assess the uncertainties of the mean 
global physical parameters and the distance, the IPM was 
run with four different inte rnal settings (for details see 
ISodor. Jurcsik fc Szeidl|[200l . Table 1) using the mean light 
curves of the two observing seasons. These runs yielded eight 
results for each physical parameter. Their averages are listed 
in Table [7] The errors of the output parameters, with the 
exception of E(B — V), were estimated as the rms of the cor- 
responding eight results. The error of the derived interstellar 
reddening depended on the calibration of our photometry, 
i.e. on the accuracy of the standard reference magnitudes 
given in Table [21 The IPM yielded an interstellar redden- 
ing value of E(B — V) — 0.19 mag for CZ Lac in accor- 
d ance with the . 28 ma g value given by the reddening map 
of ISchlegel et all (119981 ). The eight different results for the 
mass and the distance of CZ Lac were within the ranges of 
0.590-0.618 M and 1182-1205 pc, respectively. These nar- 
row ranges proved the stability of the results obtained with 
the IPM. 

Possible changes in the mean radius, luminosity and 
temperature between the two seasons were also checked. Af- 
ter fixing the constant parameters (metallicity, mass, dis- 
tance and interstellar reddening) to their values determined 
in the previous step, IPM was run again for the mean pulsa- 
tion light curves of the two seasons separately. The results 
did not show any significant variation in any of the free 
parameters between the two seasons. This finding was in ac- 
cordance with the photometric results, which showed that 
the mean BV(RI)c colours of the light curves remained the 



After determining the mass, the distance and the interstel- 
lar reddening of CZ Lac, the IPM were run for the data of 
different phases of the modulations, keeping these parame- 
ters fixed during the fitting process. In this way, the changes 
in radius, luminosity and temperature with Blazhko phase 
were obtained for both modulations in the two observing 
seasons. 

The synthetic light curves of both modulations were 
divided into 20 Blazhko-phase bins of equal width. These 
data showed practically no light-curve modulation in any 
of the bins. The light curves of the modulation phase bins 
were then fitted with the harmonics of the pulsation fre- 
quency. These fits were used as the input for the IPM. The 
results, the pulsation-phase-dependent output parameters of 
the IPM were averaged for the whole pulsation cycle, i. e. 
their arithmetic means were calculated in order to detect 
any Blazhko-phase-dependent changes. At this stage, eight 
different internal settings of the IPM were used to estimate 
the uncertainties of the results. The choices of the V ra a tem- 
plate curves and their weights during the fitting process 
were the same that we had used in our earlier analyses 
(see ISodor. Jurcsik fc Szeidll |200S| , Table 1). Furthermore, 
we used two differen t A(Vi S ,d)-V s , lnp relations (for details see 
I Jurcsik et aill2008bl ). 

The mean values of the observational data, derived from 
the photometry, and the results of the IPM are plotted in 
Fig.|S]for both modulations in both seasons. In the left-hand 
panels, the observed mean quantities are plotted as a func- 
tion of Blazhko phase: the peak-to-peak pulsation ampli- 
tude in V band, the variation of the pulsation period calcu- 
lated from the phase variation of the / p pulsation-frequency, 
and the different pulsation averages of the magnitudes and 
colours. The right-hand panels show the quantities derived 
by the IPM for each Blazhko phase: the pulsation-averaged 
mean values of the radius, the surface gravity, the abso- 
lute visual brightness averaged by magnitude and intensity 
units, the luminosity and the effective temperature. In full 
accord ance with our previous findings (|jurcsik et al.|[2008bl . 
l2009al ) , the actual values of the fixed parameters had an in- 
fluence only on the average values of those physical parame- 
ters that changed with Blazhko phase, and had no effect on 
their changes around the averages. 

As the amplitudes of both modulations decreased from 
the first season to the next, the amplitudes of the changes in 
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Figure 6. Variations of the observed mean (left-hand panels) and derived physical parameters (right-hand panels) of CZ Lac during 
the Blazhko cycles. Each of the four blocks of panels corresponds to one modulation component (/ml or / m 2) in one season. Magnitude- 
and intensity-averaged brightnesses and colours are denoted by angle and round brackets, respectively. From top to bottom, left-hand 
panels show : amplitude modulation - total pulsation amplitude; pulsation-phase or period modulation - variation of the phase of the 
/ p pulsation component (empty circles) and deviation of the instantaneous pulsation period from the mean pulsation period (continuous 
line); pulsation-averaged V magnitude; pulsation-averaged B—V colour; pulsation-averaged V—Ic colour. From top to bottom, right-hand 
panels show the pulsation averages of the following physical parameters : radius; surface gravity; absolute visual magnitude; luminosity; 
effective surface temperature. 
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most of the derived parameters decreased, too. Although the 
two modulations had rather different characteristics, apart 
from the reduced amplitudes, their behaviour remained the 
same from the first season to the second. The consistency of 
the IPM results for the two modulations in the different sea- 
sons is a strong argument for the reliability of the method. 
Taking into account that the input data of the IPM are com- 
pletely independent for the two seasons, such a consistent 
result can hardly be explained by any artefact 

In accordance with ou r previous results (|Jurcsik et al.l 
l2008bl . 120093 ; ISodorl 120091 ) . the intensity-averaged mean V 
brightness reflects the mean luminosity changes of CZ Lac 
very well. The luminosity is always the highest at the high- 
amplitude Blazhko phases. The luminosity variations con- 
nected to the first modulation (/ m i) are caused predomi- 
nantly by temperature changes in both seasons. The radius 
changes corresponding to this modulation component are 
very small and are in anti-correlation with the luminosity 
variations. At the same time, the luminosity changes con- 
nected to the second modulation (/ m 2) are induced mainly 
by radius changes in both seasons. The small temperature 
changes with this modulation are out of phase in the first 
season and are the opposite of the luminosity variations in 
the second season. 

The strength of the phase (period) modulations seems 
to have no influence on the strength of any of the detected 
variations of the parameters. 



7 DISCUSSION AND SUMMARY 

Multi-periodic modulati on has already been found in some 
Blazhko RR Lyr stars llLaCluvze et all 120041: ISodor et al 



20061 ; ISzczygiel fc Fabrvckvl 120071 ; 1 Wils. Kleidis fc Broend 
20081 ; iBenko et al.ll2010T l~ However, one of the modulations 



was always dominant, with much higher amplitude than the 
other, in these stars. CZ Lac is the first multi-periodically 
modulated RR Lyrae star with two modulations of similar 
strength. The good phase coverage of our multicolour CCD 
observations made it possible to study the nature of the 
multi-periodic modulation of CZ Lac in detail. 

The analysis of the data from two consecutive observing 
seasons revealed that the Fourier spectra of the light curves 
showed complex multiplet structures around the pulsation 
harmonics; frequencies at fc/ p ±/ m i, fc/ p ±/ m 2, fc/p±(/ m i + 

/ma), fc/p ± (/ml - ,/m2), kfp ± 2/ m l , kfp ± 2/ m2 , kfp - 3/ m2 

were detected. One of the most remarkable properties of 
the multi-periodic modulation of CZ Lac was that, similarly 
to multi-mode pulsation, combination terms of the modula- 
tions also appeared in the spectra. Modulation peaks with 
separations of half the / m 2 frequency and linear combination 
terms involving this component [kf p ± 0.5/ m 2 and kf p ± 
(/ m i — 0.5/ m 2)] were also identified. No similar frequency 
component has ever been detected in any other Blazhko star. 
No additional frequency at 0.5/ m separation has been re- 
ported even in the high-precision C0R0T and Kepler space 
missions ' data (|Chadid et alj|2010l ; iBenko et al. ||2010|). The 
appearance of the 0.5/ m 2 components makes the interpre- 
tation of the spectrum rather difficult. Taking also into ac- 
count that modulation periods of 130 d and about its half, 
67 d were detected at different epochs in one of the RRa b 
stars (V2) in the globular cluster M5 (|jurcsik et alj|2010l ). 



one might have some doubt whether / m 2 or its half is the 
true base frequency of this modulation. 

The modulation properties of CZ Lac changed signifi- 
cantly between the two observing seasons, while only very 
slight changes were detected in the mean pulsation proper- 
ties. The periods of the two modulations changed in opposite 
directions: P m i increased from 18.5 to 18.7 d, while P m 2 de- 
creased form 14.8 to 14.4 d. The amplitudes of both modula- 
tions decreased significantly from the first observing season 
to the next, and their relative strength reversed. The period 
and the V amplitude of the mean pulsation light curve de- 
creased by 4 ■ 10~ 6 d and 0.015 mag, respectively. All the 
detected changes took place within a 237-d-long interval, as 
there was no sign of any frequency or amplitude change in 
the individual seasons' data. 

The two modulation frequencies of CZ Lac were in a 
5:4 resonance ratio well within the error limit in the first 
observing season. The modulation frequencies changed to 
the second season so that they were close to the 4:3 reso- 
nance then, but the frequency ratio differed somewhat from 
the exact resonance value. The sudden transition from a 5:4 
resonance to a near 4:3 resonance between the two seasons 
hints that the two modulations prefer to be in resonance 
and their frequencies are not fully independent. However, 
the strong, opposite period changes of the components have 
made it impossible to describe the light-curve variations of 
CZ Lac uniformly with a single base-modulation frequency 
in the two seasons. There is another example of a Blazhko 
star (V14 in M5) w ith modulation periods in 4:3 resonance 
(|jurcsik et alJl20lbT ). However, the modulations with the dif- 
ferent periods were not coincidental, they were observed at 
different epochs in this star. These observations suggest that 
the detected modulation periods of Blazhko stars may not be 
unique, which makes the interpretation of the phenomenon 
even more difficult. 

Changes in the physical parameters of CZ Lac 
were determined by the In verse Photometric Method 
(|S6dor. Jurcsik fc Szei d1 l2008h . No changes in the mean pa- 
rameters were derived for the two consecutive seasons in 
accordance with the constancy of the mean colours and 
magnitudes during the whole two-years observations. Small 
changes in the mean effective temperature (T e g), mean lu- 
minosity (L) and mean radius (R) with the phases of both 
modulations were unambiguously detected. The relations 
between the physical-parameter variations were different for 
the two modulations in both seasons. Although the large- 
amplitude phases of both / m i and / m2 corresponded to 1— 
3 per cent increased mean luminosity values, these luminos- 
ity changes were induced by either temperature or radius 
changes connected to one and the other modulations. This 
pattern remained unchanged from the first season to the 
next. Taking into account the independence of the input 
data (observations of the two seasons) and their separate 
analysis, such stability of the results confirms that the de- 
tected changes are real. 

CZ Lac is an example for modulations of different char- 
acter showing up simultaneously in a single star. Conse- 
quently, it seems that the stellar structure and parameters 
do not directly determine how the mean physical parameters 
vary with the Blazhko cycle. This fact may explain why no 
generally valid relation between pulsation and modulation 
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properties of Blazhko stars has been revealed in spite of the 
efforts for establishing such a connection. 

The phase relation between the pulsation-period and 
-amplitude variations during the Blazhko cycle should be 
an impor tant parameter of the modulation according to the 
model of IStothersl (|200fj| . I2010T ) . Examining these phase re- 
lations of the stars studied so far in the Konkoly Blazhko 
Survey and the data of the two modulations of CZ Lac, we 
have not found any connection between these properties and 
any of the observed and derived (by the IPM) parameters 
of the modulation. The only regularity that have been dis- 
closed is that, in all the studied cases, the changes in the 
mean luminosity are always parallel to the amplitude varia- 
tions of the pulsation, i.e. Blazhko stars are always slightly 
more luminous at their large-amplitude phase than at their 
small-amplitude phase. 

The investigation of the long-term period changes us- 
ing O — C data showed that the pulsation period of CZ Lac 
underwent a significant decrease just around the time of our 
CCD observations in 2004 and 2005. Linear fits to the O — C 
data for the 1990-2005.5 and 2005.5-2010 intervals yielded 
pulsation periods of 0.432183 and 0.432138 d, respectively. 
The CCD observations revealed that the pulsation-period 
change was not abrupt, but lasted for several years. The 
Konkoly CCD observations were made serendipitously at 
the time when the drastic period decrease started. Our data 
showed that the modulation properties (periods and ampli- 
tudes) changed significantly in 2004-2005, when the drastic 
pulsation-period decrease of CZ Lac just began. All these 
changes were presumably triggered by structural changes 
in the stellar interior. The strong variations observed in 
the modulation properties implied that the modulation was 
more sensitive to the changes in the interior and/or followed 
these changes more quickly than the properties of the pulsa- 
tion. A similar, significant period change (an increase, how- 
ever) occurred in RR Gem jSodor et al.ll2007l . Fig. 1) around 
1938, which was also accompanied by drastic changes in the 
modulation properties. 

CZ Lac is the first Blazhko RR Lyrae star that showed 
two different modulations simultaneously with amplitudes 
of about equal strengths. If the two modulations are indeed 
independent of each other, then it contradicts any expla- 
nation of the Blazhko effect that connect the modulation 
period to the rotation of the star. 
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